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We investigated the role of non-Gaussianity (nG) plays in the enhancement of the fidelity in 
continuous-variable quantum teleportation of ideal Braunstein and Kimble (BK) protocol for coher- 
ent states, theoretically. The de-Gaussification procedure is realized through subtracting photons on 
the two-mode squeezed vacuum state (TMSVs). We find that the high fidelity always refers to a sym- 
metrical arrangement of photon subtractions on the different modes of the TMSVs. The non-Gaussian 
resources demonstrate commendable superiorities compare with the Gaussian resources only for sym- 
metrical arrangements of photon subtractions, however, the asymmetrical arrangements do not. When 
the total number of photon subtractions be a constant, we find that the optimal nG procedure prefers 
the most asymmetrical arrangement of subtractions. This characteristic is not consistent with the result 
that the highest fidelity refers to a symmetrical case. Under the same squeezing parameter, a higher 
nG might not always lead to a higher fidelity. 
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■ I. INTRODUCTION 

^ ' In recent years, quantum information and communication (QIC) technologies of continuous variables have been 
^ , discussed vastly in literatures focused on the entangled non-Gaussian optical states as communication resources [ ^ - 
^ ■ 10]. This is mainly because their strongly nonclassical properties, such as the enhanced entanglements and clear 
negativity of the Wigner functions, which are important for the improvement of the efficiency in the implementation 
of QIC [11-2 ]. In turn, to circumvent the restriction of no-go theorem [22-25] for entanglement distillation of 
\ Gaussian states, it has report that using non-Gaussian states and operations are feasible protocols [26, 27]. It has also 
been shown that the non-Gaussian states are necessary to implement universal quantum computations [28-30]. 

At present, non-Gaussian quantum states can be feasibly obtained through probabilistic operating procedures in 
laboratory [31-37]. A class of non-Gaussian quantum optical states can be realized by adding or subtracting an 
' arbitrary number of photons on the TMSVs [33-38]. However, due to the rapid decrease of the success probability, it 
[ becomes increasingly more challenging to construct large number of photon subtracted or added non-Gaussian states. 
' Recently, Namekata et al. [39] reported the experimental demonstration of a one- or two-photon subtracted squeezed 
[ state. 

. Theoretically, Dell' Anno et al. [6-9] investigated the performance of a class of two modes squeezed Bell-like en- 
tangled resources that take photon subtracted and photon added as particular instances in continuous-variable quan- 
tum teleportation of ideal and imperfect BK protocols [40]. They found that the optimized non-Gaussian entangled 
resources outperformed their counterpart Gaussian resources in the implementation of QIC. However, the amount 
^ of entanglement of the communication resources are not the only factors for optimizing the teleportation fidelity. 
H Though the same degree of entanglement, different procedures for de-Gaussification, generally, gave out different en- 
- . ." hancements of the fidelity. Yang et al. [5]. investigated the relations of the fidelity with the Einstein-Podolsky-Rosen 
correlations. And they found that in weak squeezing regions of certain resource it can carry out a high fidelity. 

In this paper, we generally investigate the performance of a class of two mode non-Gaussian entangled resources 
Eq. (1) in continuous-variable quantum teleportation of ideal BK protocol for coherent states Eq. (5). The non- 
Gaussian entangled resources are obtained by performing local photon subtracted operations a and b on the two 
modes of TMSVs, respectively. In some sense, we generalize the non-Gaussian entangled communication resources of 
the work in Refs. [ 2, ]. However, in our present work we analysis the role of nG plays in the enhancement of the 
fidelity, emphatically. Effective arrangements of photon subtraction operations for the improvement of teleportation 
fidelity are discussed. 

This work is motivated by the following reasons. Firstly, in principle, the nG is not directly related to the non- 
classical character of a quantum state. It has shown that the nG and entanglement of photon subtracted states have 
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radically different behaviors [18]. Although the fidelity of teleportation can be improved by various de-Gaussification 
procedures [2, 4-y], it is not clear what role the nG plays in teleportation protocols, and whether a higher nG always 
leads to a higher fidelity. We adopt the measure proposed by Genoni et al. [41, 42] to investigate the role of nG in 
this issue. Secondly, effective arrangements for de-Gaussification operations are considered. How the fidelity of tele- 
portation can be optimized through properly arranging several limited de-Gaussification operations? In particular, we 
analyze the asymmetrical and also the symmetrical arrangements of the operations on different modes of the TMSVs 
where the total number of photon subtractions is constant. Lastly, considering photon subtraction technologies have 
been used for various aims in experiments, and further theoretical analysis have shown that photon subtractions are 
more efficient from the energy-cost point of view for increasing entanglements of communication resources [ ], only 
photon subtraction operations are included in our present work. However, the whole processing is also suitable for 
testing photon additions and other schemes for generating non-Gaussian states using the TMSVs [5, 32]. 

This paper is organized as follows: In Section II, for self-content, we introduce the theoretical expressions of the 
photon-subtracted TMSVs, and the photon subtraction techniques in laboratories, briefly. In Section III, using the 
non-Gaussian entangled resources Eq. (1), we generally investigate the teleportation fidelity of ideal BK protocol for 
coherent states (5) . The relationships among the teleportation fidelity with the entanglement and the nG are discussed 
in Section IV. The final section. Section V, is devoted to conclusions and discussions. 



II. THE ARBITRARY NUMBER OF PHOTON-SUBTRACTED TMSVS 



Theoretically the arbitrary number of photon-subtracted TMSVs can be obtained by repeatedly operating the photon 
annihilation operators on a TMSVs. The normalized form can be written as [15, 16, 43] 

=a™&"S'(r)|00), (1) 



ted p 
ibtra( 

mini sinh^" rF/;'-™^°) (cosh 2r) (m < n) 
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where m, n are the number of arbitrary subtracted photons on the two different modes (for Alice and Bob). The 
square of the normalization constant of photon-subtracted TMSVs takes 



and Pm'^^ (x) is the Jacobi polynomial which defined as 



k=0 



k \ m ^ k 



S{r) = cxp[r(a^6^ — ah)\ is the two-mode squeezing operator with r being a real squeezing parameter, and the Bose 
annihilation and creation operators satisfy [a, ~ [b, 6^] = 1. When m = n = 0, the normalization constants reads 
l/^A^r,o,o = 1, and the states reduces to the familiar original Gaussian TMSVs [44-46] 



|TMSV) = S'(r)|00) Vl - ^ A"|nn), (3) 

n=0 

where A = tanhr. 

The photon subtraction process on TMSVs was first proposed by Opatrny et al. [ ], theoretically. In laboratories, 
photon subtraction can be done by splitting a small fraction of the squeezed vacuum for photon counting [35, 36, 47], 
which are guided into a photodetector, and by selecting the transmitted state conditioned on the detection of photons 
which can be signed as "clicks" in an on/off photodetector, e.g. the avalanche photodiodes. Although subtracting 
large number of photons leads to complex states, experiments get more challenging due to the rapid decrease of the 
success probability[18]. 



III. THE TELEPORTATION FIDELITY USING THE NON-GAUSSIAN ENTANGLED RESOURCES OBTAINED BY PERFORMING 

PHOTON SUBTRACTED OPERATIONS 

A. The teleportation fidelity 

The success probability to teleport a pure quantum state can be described through the fidelity of teleportation, 
F = tr(pin/9out), which as a measure of how close between the initial input state and the final (mixed) output quantum 
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state. In the formalism of the characteristic functions of continuous variables, the fidelity can be written as [49] 

1 



F = - d axin(a)Xout(-a), 



(4) 



where tPa = dRe{a)dlm{a). Using the Weyl displacement operator D{a) = cxp{aa'^ — a* a) , the characteristic function 
of the single-mode input coherent state |^) = D{^)\0) reads 



— — h 2ilTn{afi* 



Xin(a) ^ {fJ-\D{a)\fi) = cxp 
For the output teleported state, the characteristic function has the factorized form [50] 



Xout(a) = Xin(a)xi2(a*,a), 



(5) 
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where xi2(q!, /3) is the symmetrically ordered characteristic function of the entangled resource. For the non-Gaussian 
entangled states Eq. (1), Eq. (6) reads 



Xi2{a,l3) = {i^\D,{a)D2{f3m, 



(7) 



where Di{a) ~ cxp(aa''' — a* a) and D2{P) = exp(/3&^ — j3*h) correspond to the two modes of the entangled states, 
respectively. Using the relations of operators 
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with A" (a) = (-^)™ and A"(/3) = JL^" (af^)"' ^q- (7) can be reduced to a simple form 

xlr^(a,/3) =iVf,„,„-'e-^(l"l^+l'5|^)A™(a)A"(/3) [x(a, /3)e^(l"l^+l^l^ " 
where the characteristic function of the TMSVs Eq. (3) reads 



x(a,/3) = exp 



[ 2(1 -A2) 
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and (m, n) denotes the number of photon subtracted operations on different modes (1 and 2) of the TMSVs. According 
to the above expressions, for the arbitrary number (m, n) of photon subtractions, the fidelity for teleporting a coherent 
state based on the non-Gaussian entangled resources can be obtained, which is a function of the squeezing parameter, 
viz. F ~ F{\). The fidelity is equal for the two modes (1 and 2) of photon subtracted operations and this can be 
represented as 



(12) 



For special cases, benefit from the number expressions of the series expansions for A, we find the results which are 
listed in the A. However, because of the arbitrary order partial derivatives in Eq. (10), finding general expressions 
come across challenges. 



B. Effective arrangements of photon subtractions 

Using the results in the A, it is straightforward to discuss the effective arrangements of photon subtractions on the 
TMSVs for the improvement of teleportation fidelity. 

(i) Firstly, considering the equality of the two modes, only cases n > m (with m = 0,1,2,...) are enough. We plot 
the series in Fig. 1. It is shown that, for the same squeezing parameter, non-symmetrical arrangements lead to lower 
fidelities, while higher fidelities always prefer to sjTnmetrical arrangements of photon subtractions on the two modes 
of the TMSVs. For small squeezing parameters, non-sjnnmetrical arrangements even could not guarantee that the 
fidelities over the classical limit 1/2. A significant result is that some non-symmetrical arrangements does not exhibit 
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FIG. 1: (Color online) Fidelities as functions of the squeezing parameter with n > m and (a) m — 0; (Jo) m = 1; (c) m — 2; (d) 
m = 5. For the same squeezing parameter, higher fidelities always refer to symmetrical arrangements of photon subtractions on 
the different modes of the TMSVs. The fidelity for the Gaussian resource (m = n = 0) is attached in dashed green (dashed) lines. 



superior performances compare with the Gaussian resource (m = n = 0) which is attached as green dashed hnes (see 
Fig. 1(b) -(d) and Fig. 2). 

(ii) Further, letting the total number of photon subtractions be a constant C and from the curves in Fig. 2, it 
can be shown that high fidelities also prefer to symmetrical arrangements. We plot the series of symmetrical cases 
in Fig. 3. And it is shown that the fidelities are improved entirely, where the non-Gaussian resources demonstrate 
commendable superiorities compare with the Gaussian resources. However, fidelities are not monotonically increased 
with the number of operations. 

These properties are very similar with the characteristics of entanglement in Ref [18], where it has shown that the 
highest entanglement is achieved when these operations are equally shared between the two modes. 

IV. THE TELEPORTATION FIDELITY WTTH THE NON-GAUSSIANTTY OF RESOURCES 

In this section, we investigate the role of nG plays in teleportation protocols, where the measure proposed by Genoni 
et al. [41, 42] is adopted. 

For a general state p, the degree of nG is quantified as the difference of the von Neumann entropy between the state 
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FIG. 2: (Color online) Fidelities as functions of the squeezing parameter with the total number of photon subtractions be a constant 
C: (a) C = 2; (b) C = 3; (c) C = 4; and (d) C — 10. And the fidelity for the Gaussian resource (m = n = 0) is also attached 
(green dashed lines). 



and its reference Gaussian state 

5[p\ = S{p\ \t) = Tr[plnp] - Tr[plnr] = S{t) ~ S{p). (13) 

The reference Gaussian state t has the same first- and second-order moments of canonical quadrature variables with 
the state p. The two moments are the vector of mean values X[t] ~ X[p\ and the covariance matrix (CM) o-[t] = (^[p], 
respectively, which defined as 

X := {R), (14) 
+ (15) 

where the vector operator i? := ■ • • , x„,p„)^. For the two-mode states, the formula can be reduced to 

5[p] = h{d^) + h[d+)-S{p), (16) 

and the function h{x) is defined as 

h{x) = {x + \)\n{x + \)\n{x - 1). (17) 



FIG. 3: (Color online) Fidelities as functions of the squeezing parameter for symmetrical arrangements of photon subtractions. 
Although fidelities are not increased monotonically, they are being improved entirely compare with the Gaussian resources (m = 
n = 0). 

Here d± are the symplectic eigenvalues of the two-mode CM 



A 


C 




B , 



(18) 

with A, B, and C are 2x2 real matrices. In the standard from, matrices A and B are proportional to the identity and 
C is diagonal. Using the four local sjnnplectic invariants Ii = det(A), I2 = det(i?), /a = dct(C), and = dct(o-), the 
symplectic eigenvalues can be expressed as 



C^-(n+l)sinh(2r) - „^ . , M . J- (22) 



. ^m±^^E^ (19) 

with d+ > c?_ and A{a) ~ Ii + I2 + 2/3. Actually for pure stated p since S{p) = 0, the nG is always equal with the 
von Neumann entropy of the corresponding Gaussian states t. Using the definition of CM Eq. (15), one can obtain 

^=(^ + %^)a, (20) 

S=(^ + %^)a, (21) 

, .^....^,.o„^ -P"""'''(cosh2r) /l 
Pi"-™'"'(cosh2r)U -1 

Without loss generality, we have assumed that 71 > to in the above results. Now, it is straightforward to give the nG 
for pure states Eq. (1) according with the definition in Eq. (16). 

In Fig. 4(a)-(c), we plot the nG as the functions of the squeezing parameter. For cases with m ■ n ^ 0, the nG is a 
monotonically increasing function. And for series with n — m = Constant, they have the same nG as the squeezing 
parameter comes to zero, see Fig. 4(a). The nG for series n + to = 10 are ploted in Fig. 4(b). It is shown that the 
more asymmetrical arrangements of the operations on the two modes always prefer to higher nG, which is consistent 
with the findings in the literature [18]. When photon subtracted operations are only added on the one mode (e.g. the 
most asymmetrical arrangement to = 0), the nG reaches a maximum S[p] = ~n\n{n) + {n + l)ln(77, + 1) (it becomes 
does not rely on the squeezing parameter) . However, the optimal procedure for getting a highest nG prefers the most 
asymmetrical arrangement for the operations when n + m ~ Constant, is not consistent with the result that the highest 
fidelity always prefers to a symmetrical arrangement in Sec. IIIB. For the series of symmetrical cases, it can be seen 
in Fig. 4(c) that a larger number operations always represent a higher nG. In Fig. 4(d), we plot the relations of the 
fidelity with the nG for sjTnmetrical cases. Though all the fidelities are being raised due to the introduction of the nG, 
the case of to = n = 1 represents the most rapidly increasing behavior 
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FIG. 4: (Color online) The nG as the functions of the squeezing parameter in figures (a)-(c). (a) When n — m = Constant, the 
series have the same nG as the squeezing parameter comes to zero, (b) More asymmetrical arrangements of the operations on the 
two modes give more nG. (c) For symmetrical arrangements, a larger number operations always represent a higher nG. Figure (d) 
shows the relations of the fidelity with the nG for symmetrical cases, where m = n = 1 represents the most rapidly increasing 
behavior. 



V. CONCLUSIONS AND DISCUSSIONS 

In this paper, we studied the effects on the fidelity of a de-Gaussification procedure on the entangled Gaussian com- 
munication resources in continuous-variable quantum teleportation of BK protocol for coherent states. The procedure 
is realized through subtracting photons on the TMSVs. We investigated the role of nG plays in the enhancement of 
the fidelity, emphatically. 

We find that the high fidelity always prefers to a symmetrical arrangement of photon subtractions on the different 
modes of the TMSVs. When the total number of photon subtractions is a constant, the conclusion is still affirmative. 
Finding a general expression of the fidelity as a function of the squeezing parameter comes across challenges due to 
the arbitrary order partial derivatives. Using the mathematical induction method, we find the analytical results for 
some special cases which are listed in the A. 

A significant result is that non-Gaussian resources demonstrate commendable superiorities compare with the Gaus- 
sian resources only for symmetrical arrangements of photon subtractions. Non-symmetrical arrangements even could 
not guarantee that the fidelities over the classical limit 1/2. They do not exhibit superior performances compare with 
the Gaussian resource. 
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The role of the nG plays in teleportation protocols is studied. When operations take n + m = Constant, the optimal 
de-Gaussification procedure prefers the most asymmetrical arrangement (m = 0) . This characteristic is not consistent 
with the result that the highest fidelity prefers to a symmetrical arrangement. In symmetrical cases, a larger number of 
operation represents a higher nG, however, the lowest case m = n = 1 exhibits the most rapidly increasing behavior 
Then at the same squeezing parameter, a higher nG might not always lead to a higher fidelity. 

Considering the recent study in Ref. [18], where the authors invested the quantum entanglement properties en- 
hanced by photon addition and subtraction, in our work we find that the characteristics of the improved fidelities are 
more preferred to the enhanced quantum entanglement. While, the relations between the improved fidelities and the 
nG, which measured by the method of Genoni et al. [41, 42], are subtle. In certain conditions, enhancing the nG is 
not always means improving the fidelity of the quantum teleportation. 

The weakness of the present work is that we do not consider the rate of success of the probabilistic photon subtrac- 
tion operations performing on the initial entangled TMSVs. In the very recent and interesting work [ ], the authors 
analyzed different strategies of local photon subtraction from TMSVs in terms of entanglement gain and success prob- 
ability. For the lossless limit case, where the loss of the fraction of photons is zero, they found the best strategy to 
produce the highest entanglement gain rate is subtracting a single photon only from one mode. Affirmatively, from 
our asymmetric cases in Fig. 1(a), one can find this strategy will do not work well since it even could not guarantee 
the fidelities over the classical limit 1/2. For the lossy case, they found the best strategy is symmetric subtractions. 
Interestingly and fortunately, our result also shows that the high fidelity always prefers to a symmetrical arrangement 
of photon subtractions. However, a thoroughly analysis is interesting and needed, which is beyond the present work 
and will be discussed further in the future research directions. 

Acknowledgments The author thanks Ph.D San-Min Ke for very helpful discussions. This work was supported by 
the Basic Research Foundation of Engineering University of CAPF(WJY-201104),and partly supported by the National 
Natural Science Foundation of China under Grant Nos. 61072034. 

Appendix A 

Due to the arbitrary order partial derivatives in Eq. (10), finding general expressions of the results of Eq. (4) comes 
across challenges. For some special cases, benefit from the number expressions of the series expansions for A, we find 
the results can be written as the forms which are listed in the the following. 

Defining the function as 



y(m,n) _ 



A2" m!7i!(l + A) 

AjS ')m+n+l(1 _ \\m+n^ 



(Al) 



we have 

(i) For m = and n>m, the fidelity takes 



(A2) 



For resouces of the TMSVs, it is simpliy reduced to F{X) = /(°^"' = (1 + A)/2. 
(ii) For m = 1 and n > m, the fidelity takes 



(iii) For m = 2 and n> m, the fidelity takes 



-n ~ 4nA + {n + 1)A^ 



- 1) - yn(n - 1)A + 47i(3?i - 1)X^ - 4n{n + 1)X^ + l!l±llj!i±^ 



f(2,n) 



(A3) 



(A4) 



(iv) For m = 3 and n>m, the fidelity takes 



p{3,n) 



26 



— 7i(7i - l)(7i - 3) -n{n - l)(7i - 3)A + 8n{n - l){5n - 7)X'^ - —n{n - l){5n - 1) 



3! 



+ 2n{n + l){bn - 2) - 2n{n + l)(n + 2)A^ + 
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(v) For m = 4 and n > m, the fidelity takes 
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(vi) For m = 5 and n>m, the fidehty takes 
"210 
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^1 n(n - l)(n - 2)(n - 3)(n - 4) - -^n{n - l)(n - 2)(n - 3)(n - 4)A 
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+ - l)(n - 2) (62 - 77?i + 2171^) A'' - ^n{n ~ \){n - 2) (6371^ - 9l7i + 6) A^ 

4 8 
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O O 
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